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I discuss how a variatonal approach can be extended to systems of identical parti- 
1^ ■ clcs (in particular fermions) within the path-integral treatment. The applicability 

of the many-body variational principle for path integrals is illustrated for different 
model systems, and is shown to crucially depend on whether or not a model system 
tt , possesses the proper symmetry with respect to permutations of identical particles. 
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1 Introduction 



In the path-integral formulation of quantum mechanics, the so-called Jensen- 
Feynman inequality provides an upper bound tCLthe free energy of a quantum 
system, if properly applied. It was introduceda in Feynmaii's path-integral 
approach to the Frohlich polaron (see formula (8.40) in Ref.a): 

F < F M + ~<S - S M )s M if S, S M are real. (1) 

\o ■ p 

In the variational functional, F and S are the free energy and the action 
functional of the system under consideration, whereas F M and Sm are the 
i-v , free energy and the action functional of a model system; the temperature is 

described by the parameter (3 = l/fceT. Angular brackets mean a weighted 
average over the pathsa: 

ux „ _ J»exp(-5 M )£>(path) 



/ exp(-S , M )-D(path) 



G 

A rigorous argument to prove the inequality (|lj) is based on the convex nature 
£h I of the. exponential exp(x) of a real stochastic variable x (see e.g. Fig. 11-1 in 

Ref. H), which leads to (e x ) > e^with (x) the weighted average of x. 

Apart from Feynman's variational treatment of the ground state energy of 
a polaron, the path-integral approach based on the Jensen- Feynman inequality 



a It is implicitly assumed that the action functional and the path integral are expressed in 
the imaginary-time variable. This convention is followed throughout the present paper. 



was successfully applied to a series of .CEpblemsO, e.g., to the calculation of 
the effective classical partition fjmctioroEl and of quantum corrections to the 
free energy of nonlinear systemsEm, to the description of all critical exponents 
observahlfipLo second-order phase transitions^, and to the problem of bipolaron 
stabilityBB. 

The derivation of the Jensen-Fcynman inequality crucially depends on 
the assumption that both the action and the trial action are real functionals. 
As already recognized by Feynman (see Ref.tll, p. 308), its application to a 
polaron in a magnetic field therefore becomes problematic, because the action 
functional S for a polaron in a magnetic field (and any reasonable trial action 
Sm) is no longer real- valued. A discussion of this problem lies beyond the 
scope of the prpspnt .paper. Fy rnnrp details on the status of this problem, 



pe ot tne present-pa 
see the literatureBEaM 1151 ' ' : 



The Jensen-Feynman inequality is reminiscent of the Bogolubov 
inequalityEjtJ, which provides the following upper bound to the free energy 
F of a system described by the Hamiltonian H 

1 Tr \(H - H M ) e -P HM ] 
F < F M - - lK M > 1 if H, H M are Hermitian, 3 

p fr(e _f " lM ) 

where Hm is the Hamiltonian of some trial system_with free energy Fm . The 
Rayleigh-Ritz variational principle (see e.g. Ref£d, p. 172) for the ground 
state energy E < (\P M \H \ ^> M ) / (*m|*m) with a trial state |* M ) is the 
zero-temperature limit of the Bogolubov inequality. 

The condition that S and Sm are real in (fil) is not necessarily equivalent to 
the requirement that H and Hm are Hermitian in (ph . If the Hamiltonians H 
and Hm in the Bogolubov inequality are Hermitian operators corresponding to 
Lagrangians L and Lm in the Jensen- Feynman inequality, then the one-to-one 
correspondence between ffl) and (H) guarantees the validity of the Jensen- 
Feynman inequality, even if the action functionals are not real (e.g. for a 
particle in a magnetic field). 

However, both inequalities do not necessarily have the same physical con- 
tent: for a system with action S it is not always possible to derive a corre- 
sponding Hamiltonian. For example, for the Frohlich polaron (in the absence 
of a magnetic field) the strength of the Jensen-Feynman inequality lies in 
the fact that it remains valid after the elimination of the phonons, with a 
retarded effective action functional, for which no corresponding Hamiltonian 
representation is known. In the operator formulation, the phonon elimination 
can formally be realized with ordered-operator calculus, but this approach 
involves non-Hermitian effective operators in the electron variables. 

Fermion systems (with parallel spins) form an important class of systems 



for which the Jensen-Feynman inequality is not directly applicable (whereas 
the Bogolubov inequality remains valid in the Hilbert space of antisymmetric 
states under permutations of the particle coordinates). The reason is that the 
path integral for fermions with parallel spin, if expressed in the full coordinate 
space, is a superposition of path integrals with all possible permutations of 
the particle coordinates, with negative signs for all odd permutations. For 
bosons, no negative signs result from the permutations and the application 
of the Jensen-Feynman inequality presents no problems. Therefore, only the 
many-fermion problem will be explicitly addressed below. 

2 Path-integral approach for many-body systems 

Recent studies on the path-integral approach to the many-body problem for 
a fixed number of identical particles by Brosens, Lemmens and Devreesec3 
have allowed to calculate the Feynman-Kac functional on a state space for 
N indistinguishable particles, which was found by imposing an ordering on 
the configuration space, and the introduction of a set of boundary conditions 
in this state space. The path integral (in the imaginary-time variable) for 
identical particles was shown to be positive within this state space. This 



implies (see subsection 3.1 for more details) that a many-body extension of 
the Jensen-Feynman inequality was found, which can be used to evaluate The 
partition function for interacting identical particles (Ref.E3, p. 4476, Ref.El). 
This many-body variational principle for path integrals was applied to the 
study of thermodynamic^], properties of a spin-polarized gas of bosons (Refo, 
abstract, Eq. (3); Ref.c3, Eq r -.(13)). The applicability of the variational 
principle as formulated in RcfE3 for manypbody problems, was discussed in 
relation to the analysis of correlations (Ref.Ea, p. 1641) and thermodynamical 
properties (Ref.EZl, p. 3911) of a confined gas of harmonically interacting 
spin-polarized fermions. 

The many-body variational principle for path integrals was also used re- 
cently in order to calculate the ground state energy and the optical absorption 
spectrum of a many-polaron system, confined to a quantum dot (Ref.c3, p. 
306). 

The remainder of this paper addresses the question which choice of model 
actions is allowed in order to treat specific systems of interacting bosons and 
fermions. I will give some examples, illustrating that the applicability of 
the many-body variational principle for path integrals crucially depends on 
whether or not a model system possesses the proper symmetry properties with 
respect to permutations of identical particles. The requirements analyzed in 
this article are qualitatively new as compared to the Feynman variational 



principle of Refs.fl 

3 Many-body variational principle for path integrals 

Let a many-fermion system be described by the action functional 5[x(£)], 
where x = {xi, . . . , x^v} are the coordinate vectors of fermions. The partition 
function Zp of a many-fermion system can be expressed as a path integral: 

^ = E^p/ dS j ( * US (t)exp{-S[X («)]}, (4) 

where the summation is over all elements P of the permutation group. The 
weight (—1) is the character of the representation, i.e. +1 for even permu- 
tations and —1 for odd permutations (for the case of fermions). 



3.1 Model systems with local potentials 

In Ref.E3, a many-body problem was analyzed for a local potential V(x) (in- 
cluding interparticle interactions) with the action functional in the imaginary 
time representation 



1 [** 

S[x(i)] = -/ dt 



N 

2£>?(t)+v(«(t)) 



(5) 



Note that this action is invariant under the permutations of any two fermions 
at any (imaginary) time, since the potential can not make a distinction be- 
tween identical particles. 

If the potential V(x(t)) is invariant with respect _to_the permutations 
of the carthesian components of the particle coordinatesE3'cB, the many-body 
propagator was obtained by four independent processes^er pair of particles, 
defined on a state space (D^ in the notations oi-.Ref.E5) with well-defined 
boundary conditions (see Eqs (4.16-4.17) in Ref.c3 for details). These pro- 
cesses were shown to give positive contributions to the propagator. Hence, 
the propagator itself is positive on D„, implying that the Jensen- Feynman 
inequality can be used to estimate the partition function for interacting iden- 
tical particles[} The partition function can then be represented (apart from a 



6 If the potential is only invariant under permutations of the particle coordinates, the sub- 
processes are not linearly independent, and transitions between the subprocesses have to 
be taken into account. The actual analysis in terms of the state space D^ is then only 
feasible in practice for a very limited number of fermions. In this case, an overcomplete 
space covered by all even permutations of the particle coordinates is more appropriate. 



normalizing factor) in the form 

Z F = [ dx [ D5L{t)exp{-S[5i{t)]}, x(t) £ D\ 
Jdi Jx 



(6) 



Consider now a model system with the action functional in the imaginary 
time representation 



S M [*(*)] 



1 



h S 



(It 



N 



£x|(*)+V M (*(*)) 



(7) 



where the model potential Vm (x) contains some variational parameters, and 
allows for an analytical calculation of the path integral. Suppose furthermore 
that it is invariant with respect to the permutations of the components of the 
particles positions. The path-integral expression for the partition function of 
the model system is thus: 

Z M = [ dx / Dx(i)exp{-S A f [*(*)]}, 5t(t)eD*. (8) 

JDI JR 

One can represent (g) as follows: 

Z F = f dx f Dx (t) cxp {-S M [x (t)] - {S [x (t)] - S M [x (t)])} = 

= Z M (exp{-(5[x(*)] - 5m [x(t)])}) SM . (9) 

Here, the angular brackets denote the quantum statistical expectation value: 

(•)<?„ = [2m]" 1 / dx / £>x (t) . exp {-5 M [x (t)]} , (10) 

analogously to Eq. 0). The key element of this definition is that the path 
integrals in (Jq-|l0|) are defined on the same state space D^ 7 which stems from 
the symmetry properties of the true action S [x (£)] . 

Taking into account that the propagators are positive on the domain Df. t , 
one obtains the inequality 

Z F >Z M exjp{-(S[5t(t)}-S M [x(t)])s M }, (11) 

which is readily converted into an upper bound for the free energy 



F F <F M + ^(S[5t(t)}-S M [x(t)])s M . 



(12) 



This many-body variational principle for path integrals is formally very sim- 
ilar to the Jensen- Feynman inequality (yj). The difference between Eqs. ( |l2| ) 



and (p lies in the definition of the expectation values. In (jll]) and (|1J) the 
expectation value (|o|) is defined over a subdomain D^ of the configuration 
space, whereas the expectation value (0) in (pi) is defined over the full con- 
figuration space. However, because the symmetry properties allow to unfold 
the state space into the full configuration space, the restriction to the state 
space can be omitted in the calculation. The state space (D^ in this example) 
only serves the goal to check whether the action and the trial action have the 
correct symmetry properties. 

3.2 Model systems with retarded effective interactions 

The action functional S[x(t)] of the system under study can contain a re- 
tarded effective interaction. This is the case, e.g., for a system of N polarons 
after the phonon variables have been integrated out. Such many-fcrmion sys- 



tems substantially differ from those considered above in subsection 3.1, and 
a different class of model systems seems appropriate. For this purpose, we 
consider a model system consisting of fermions in interaction with auxiliary 
fictitious particles. 

Such a model system has the action functional in the imaginary-time 
representation 

1 f hf3 
Sm = j. L M (t)dt. (13) 

The model "Lagrangian" is chosen in the form: 

L M = L F (x) + L f (y) + L F _ f (x,y), (14) 

where {x^} = x are the coordinate vectors of the fermions, and {yj} = y are 
the coordinate vectors of the fictitious particles. The "Lagrangians" Lp(5c), 
Lf{y) and I/i?_/(x, y) describe fermions, fictitious particles and the interac- 
tion between the fermions and the fictitious particles, respectively. Here, the 
discussion is limited to the case of distinguishable fictitious particles for the 
sake of simplicity. 

The partition function Zm of the model system can be written as the 
following path integral: 



z .w = Yl |n / M J d? / m (t) / Df (t) cxp { ~ Sm} ■ (15) 



p 



Integrating out the coordinates of the fictitious particles, the partition func- 



tion (|lj) takes the form 

Z M = Z Z f , (16) 

Z ^ = zZjfiJ d ^J_ ^x(t)exp{-5 [x(i)]}, (17) 

where Zf is the partition function of the system of fictitious particles: 

Zf = JdyJ*Dy(t)expl~J o L f (y)\, (18) 

and So [x (£)] is an effective action which only depends on the fermion variables 

So[x(t)} = -J (L F (t)dt + <S> [5L(t)}). (19) 

The last term in ( pT9| ) is referred to as the influence phase of the fictitious 
particles. It is defined as 

cxp{-$ [x (£)]} = 

= [Z f ]- l Jdyj y Df(t)expl~J o [L f (y) + L e _ / (x,y)] dt I . (20) 

For interaction Lagrangians Lf_j(x, y) which are quadratic in y, the influence 
phase can be shown to take the form of a retarded effective interaction: 

ft/3 ft/3 

$ [x(t)] = I ' dt I dsK(t, s)X (t) • X (s) , (21) 



where K(t,s) depends on two time variables (t,s), while X (t) is a linear 
function of the fermion coordinates x (£) (see the next section for specific 
examples) . 

If the action functional S [x (i)] of the system under study satisfies the 
permutation symmetry conditions discussed in the previous subsection, its 
partition function can be represented as a path integral over the space state 
D\ - in the form ^ : 

Z F = f dx / Dx(t)exp{-5o[x(t)]-(5[x(t)]-5o[x(t)])}, 

= Z (exp{-(S[5t(t)}-S [5t(t)})})so- (22) 



If the model action So [x (£)] also possesses the above symmetry properties 
with respect to permutations, its partition function ( |l7| ) can also be written 
in the form of a path integral over the domain _D^, 

Z = [ rfx / Dx(i)exp{-5o[x(t)]}, (23) 

JD\ Jx 

and a quantum statistical expectation value can be defined in (123) 

(.) 5o = [Zo]- 1 f dx f Z»x (t) • exp {-S [x (*)]} , (24) 

analogously to the definition (|2|). The fact that the propagators are positive 
in the integration domain £>„ guarantees that the inequality 

Z F >Z cxp{(-(S[5i(t)}-S [5i(t)})) So } (25) 

holds true. Consequently we obtain an upper bound for the free energy, similar 
to Eq. ©: ' 

F F < F v = F + ~<S [x (*)] - S [x (t)]) So . (26) 

where Fq = — 4 In Zo and -FV = —4 In -Zp. 

Because of the presence of the retarded action (gl]) resulting from the 
elimination of the fictitious particles, one should guarantee that the functional 
So [x (£)] has the required symmetry with respect to permutations which allows 
that the many-body processes, related to the quantum statistical expectation 
value (pj), are restricted to the state space D^ at any time. This condition is 
an essential ingredient for the justification of the many-body variational prin- 
ciple for path integrals (|13). Like in the case ( |12| ) of local potentials, the upper 
bound (J2fl) to the free energy is formally very similar to the Jensen-Feynman 
inequality ( p]): the difference lies again in the definition of the expectation 
values. In (E3) and (Eff) the expectation value (p3) is defined over the sub- 
domain D„ of the configuration space, whereas the expectation value (0) in 
(uh is defined over the full configuration space. However, like in the previous 
subsection, the state space (D^ in this case) is only needed to check whether 
the symmetry of the action and the trial action allows to apply the inequal- 
ity. The calculation can be performed over the total configuration space by 
unfolding the state space. 



4 Examples: non-interacting fermions as a test case for a 
many-body variational principle with path integrals 

4-1 Model system, in which each fermion harmonically interacts with 
one fictitious particle 

In order to illustrate the applicability of the many-body variational principle 
for path integrals (|26|), and in particular the need of the correct symmetry 
requirements for the model action, we first consider a very simple system of 
N = X)cr=±i/2 N a non-interacting fermions, described by the Lagrangian 

^ = y E E*L> (27) 

er=±l/2j'=l 

where N a is the number of electrons with spin component a = ±1/2. The 
ground state energy for the system with classical "Lagrangian" ( p7| ) is ele- 
mentary: 

with the Fermi wave number kp and the Fermi energy Ep. 

We now examine whether the many-body variational principle for path 
integrals (Eq) indeed provides an upper bound to the correct ground state 
energy ([28|). 

For the model system (O) we choose a Lagrangian Lm, in which each 
fermion harmonically interacts with one fictitious particle. We do this uncrit- 
ically, deliberately overlooking the problem of the required symmetry of the 
state space of the model action and choose 

_ N " 

Lm = y e e*l+ 

o=±l/2 j=l 

+ t E E*k+£ E Efe-y.) 2 . CM) 

cr=±l/2j=l cr=±l/2i=l 

where M is the mass of the fictitious particles, and k is the force constant of 
the elastic bond between a fermion and its accompanying fictitious particle. 
We introduce the following notations: 

[k mM 

v = \ / — , u = . 

V u m + M 




For this particular case, the elimination of the fictitious particles leads to the 
influence phase (EG): 



*o [*(*)] = - 



Mw 3 
8,h 



h/3 h/3 



cosh-i 



dt I ds- 



(I 



M 

2 



sinh i hf3' 



w 



tr=±lj=l 



A* to] 2 , 

(30) 



with a quadratic effective retarded self interaction for each fermion. 

Then we apply the many-body variational principle for path integrals ( p6| ) 
naively for the chosen model system. The free energy F v from this inequality is 
calculated analytically. The parameters M and k of the model "Lagrangian" 
( |l4| ) are then found by minimizing the value of the supposed upper bound F v . 
This calculation has shown that the many-body variational principle for path 
integrals ( f26[ ) is violated for this model system, as clearly illustrated in Fig. H. 

What was wrong in the above approach? The answer is immediate: "Of 
course, we forgot to check whether the model action has the required symme- 
try!" The model "Lagrangian" ( |29| ) is not symmetrical with respect to the 
permutations of the fermion coordinates Xj, CT , because each of them is linked 
with a particular fictitious particle. 




Figure 1. _Eree energy F v (0) (6 = w/v) for non-interacting fermions [according to the in- 
equality (Eg)] (solid line) compared with the exact free energy F B (dashed line) at the 
optimal value of the variational frequency parameter v = 1.83 (in units of Ep/h, where 



Ep = 49 meV is the Fermi energy for the density n e 



5 X 10 19 cm 



chosen); = 16.4 (in units of E F ). The mass of a fermion is taken m 
electron mass. 



rather arbitrarily 



'«o 



the bare 
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4-2 Model system, in which each fermion has equal elastic bonds with all 
fictitious particles 

As a second example, we study the many-body variational principle for path 
integrals (Eft) for N — ^ J _ ±1 , 2 N a non-interacting electrons in a parabolic 
confinement, described by the Lagrangian (in the imaginary-time variable) 



L 



IT) 



N a 



E Eft 



cr=±l/2i=l 



(31) 



In particular, the case fio — > will be analyzed (a translationally invariant 
system, equivalent to the free-fermion model in subsection 4.1). The exact 
ground-state energy of this system can immediately be written down: 



tr=±l/2 ^ 



«o (N a ) + - 



^N Q {n (N a ))[n {N a )+A] 



(32) 



Here, no(N a ) is the number of the upper fully occupied energy level for N a 
fermions with spin component a . The number of fcrmions in all closed shells 

[JVoMiVc)) < N a ] IS 



A^o (no) = £ 



(n + l)(n + 2) I , s , w N , s 

^ ^ ^ = -(no + l)(no + 2)(n + 3). (33) 



A model system is now considered, which consists of particles in a har- 
monic confinement potential with elastic interparticle interactions as studied 
in Ref.o. The "Lagrangian" of this model system is chosen in the form 



N* 



N„ 



N„ 



^ = y E E(*L+^ 2x 

<7=±l/2j=l 



2 ) 



E E E E( x j> _x ^) s 



+ YE(y|+^l) + ^EE(y,--^) 

z=i i=i i=i 



o-=±l/2i=l t=±1/2 (=1 

2 



9 E EE^-y) 2 



(34) 



<r=±l/2i=l «=1 



The frequencies $7, w, f2s, w_b, the mass M of a fictitious particles, and the 
force constant fc are treated as variational parameters. It is important to 
stress, that in this model system each fermion has identical elastic bonds with 
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all fictitious particles, and therefore permutations of the fermion coordinates 
leave the "Lagrangian" ([34|) invariant. 

After integration over the paths of the fictitious particles, the partition 
function (Rq) becomes 



Z f = 



sinh 



npn L 



sinh 



/ hfllVB \ 



with n B = \ n 2 



kN F 



1 3 



w B = \ W-N B U) 2 



-\3N B -3 



V 2 



" 



kN f 



m V m 

The action (|19|) in the imaginary time representation takes the form 



(35) 
(36) 



S 



F0[X 



(*)] = 



with X 



Sro[X(t)] + *o [*(*)], 


(37) 


U dt 


^ 2 at2 

v-—n v~^ in .n , . muz I\ 9 

E E 2 *U*) + 2 x 


(38) 




o-=±l/2j=l 




N„ 




E E x ^< 


(39) 


tr=±l/2j=l 







and the "influence phase" (BOJ) for the present model becomes: 

^b ( \t — S 



$o [x(i)] 



ft/3 ft/3 -, 

fc 2 7V 2 Ar£ /-,/■, cosh 
at / as 



" B 



4m B hQ B 



M. 

2 











sinh 



I hjKiB \ 



X(i)-X( S ). (40) 



As distinct from (E9), this model Lagrangian (£34) is invariant with respect 
to permutations of the components of all fermion coordinates. Hence, for the 
chosen model system the sym metry conditions on the action are fulfilled (as 
formulated in subsection 4.1), which ensures the validity of the many-body 
variational principle for path integrals (|26|). 

The functional [resulting from ( P6[) at zero temperature] for the ground 
state energy of N fermions in a parabolic confinement with the confinement 
frequency f2g takes the form 



E v (fli,fl 2 ,w,w B ) = n 



ni 



2w 2 



E(w,N) w 



■(fii 



2 2 

at 



±(nZ-nl-nl + w B )J2%: 



37 2 



E 



f-f f^ 4w B f-f fij (fij + KJfl) I 

2—1 2—1 V ' ) 



&>2 - Wb) 

(41) 
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where the following notations are used: 



til-wi\ 1/2 /lui-n^ 1 -' 



ai =Ur4J ' a ^\w^i) ' (42) 

7 =[(^- W |)K-^)] 1/2 , (43) 

E{w,N) = 

= hw £ {iV CT [no (iV CT ) + §] -|AT (no(iV <r ))[no(iV (7 ) + 4]}. (44) 

o=±l/2 

The difference between the upper bound to the ground-state energy (Ell) 
and the exact ground-state energy (B2h is clearly positive: 



with 



E v (n 1 ,n 2 ,w,w B )-E (fl ,N) = ^ n ° ^ (45) 



and it follows from (|4^) that the inequality 

E v (n 1 ,n 2 ,w,w B )>Eo(no,N) (46) 

holds true for any values of the variational parameters, in accordance with the 
many-body variational principle for path integrals (pq). 

The minimal value of the functional (J4l|) , which is achieved at z = Qqi co- 



incides with the exact ground-state energy (B2j) of N non-interacting fermions 
in a parabolic confinement with the frequency Slo- This result confirms the 
applicability of the many-body variational principle for path integrals (Eq) to 
the many-fermion system under consideration, with a model system, whose 
"Lagrangian" (B41) is symmetric under the permutations of the components of 
the fcrmion positions. 

Thus, the many-body variational principle for path integrals ( p6[ ) is sat- 
isfied for a system of non-interacting fermions in a parabolic confinement po- 
tential (including the translationally invariant case Qq = 0), when the model 
system with the "Lagrangian" (M) is considered. 

Clearly, the Lagrangian (|31|) was not considered for its own sake, since 
it is trivial to treat. It was merely presented as a test case to illustrate our 
many-body variational principle for path integrals with two non-trivial trial 
actions which elucidate the crucial role of the correct symmetry requirements. 
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5 Conclusions 

Summarizing, the applicability of the many-body variational principle for path 
integrals (Eq) crucially depends on the symmetry of the model system with 
respect to permutations of identical particles. 

The invariance of the action functional 5o[x(t)] and of the model ac- 
tion with respect to permutations of components of the positions of any two 
fermions at any (imaginary) time ensures that the fermion propagator is pos- 
itive on the state space D^- 

It should be noted that these rather stringent symmetry conditions are 
used as an example. A more general variational principle for identical par- 
ticles, not limited to D^ but to a much larger subspace of the configuration 
space, will be presented in future publications. 

The main result of the present analysis is that a many-body variational 
principle for path integrals ( |26|) can be found in the framework of the many- 
body path integral approach, even for retarded effective interactions, provided 
that the model action and the true action have the appropriate symmetry 
properties under permutations of the particle coordinates. 
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